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1 Introduction and analysis method 

Differences between the properties of gluon and quark jets have been convincingly 
established by experiments operating at the e + e _ collider LEP at CERNlU. The most 
conclusive results have been from so-called Y events which are three jet events in 
which the angle between the highest energy jet and each of the two lower energy 
jets is about 150°. Table [l] summarizes the results for the mean charged particle 
multiplicity ratio between gluon and quark jets, r ch , found using Y events at LEP. A 
60% variation in the value of r c h is observed if the Cone or JADE-E0 jet finders are 
used to identify the three jet events instead of the k± (“Durham”) jet finder. 

QCD analytic predictions, .exist for the ratio of the mean particle multiplicity 
between gluon and quark jetscHa. A quantitative test of these predictions has not been 
possible, however, because of differences between the theoretical and experimental 
definitions of the jets and event samples. The analytic calculations employ definitions 
of the event samples and jets which are entirely inclusive. For the calculations, two 
samples of events are chosen: a sample of gluon-gluon gg events produced from a color 
singlet point source is used to define the gluon jet properties and a sample of quark- 
antiquark qq events produced under the same circumstances is used to define the 
quark jet properties. The gluon and quark jet characteristics are given by inclusive 
sums over the particles in these two samples. Thus, the theoretical results are not 
restricted to three-jet events and do not employ a jet finder to assign particles to the 
jets, in contrast to the Y event studies and other previous experimental studies of 
gluon and quark jet differences at high energy colliders. 

The experimental difficulty in obtaining a jet definition corresponding to the the¬ 
oretical one lies in the gluon jet sample, since gg production from a point source does 
not occur naturally in e + e _ annihilations. In contrast, the qq sample employed by 
the theory is the inclusive e + e‘ multihadronic one and so is readily available. In ref. 0, 
a method was proposed for LEP experiments to identify gluon jets using an inclusive 
definition similar to that used for the analytic calculations. The method is based on 
rare events of the type e + e - ^ qqginci. in which the q and q are identified quark jets 
which appear in the same hemisphere of an event. The quantity gi nc i., taken to be the 
gluon jet, is defined by the sum of all particles observed in the hemisphere opposite to 
that containing the q and q. In the limit that the q and q are collinear, the gluon jet 
ginci. is produced under the same conditions as the gluon jets in gg events from a color 
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Experiment 

k± Cone JADE-EO 


ALEPHU 

DELPHI0 

opal! 

1.19 ±0.04 

1.24 ±0.03 — 1.37 ±0.04 

1.25 ±0.04 1.10 ±0.03 


Table 1: LEP results for the charged particle multiplicity ratio between gluon and quark jets, r c h., 
for Y events defined using the k ±, Cone and JADE-EO jet finders. 


singlet point source. The jets gmci. therefore correspond closely to single gluon jets 
in gg events, defined by dividing the gg events in half using the plane perpendicular 
to the principal event axis. 

Although gg events from a color singlet point source are not produced in e + e - 
annihilations, they may be generated using QCD Monte Carlo programs such as the 
Herwig £3 and Jetset liTl parton shower models. To illustrate that the definition of 
gluon jets gi nc i. from the e + e _ events corresponds closely to that employed by theory 
from the gg events, the Monte Carlo is used to calculate the mean charged particle 
multiplicity in gg events, (n c h.) gg events, and the results are divided by two so that 
they correspond to a single hemisphere: these results are shown as a function of 
the jet energy Ej et =E c . m ./2 by the solid and dashed lines in Fig. [l](a). Shown in 
comparison, by the cross and diamond symbols, are the results for g inc i. jets from 
e + e _ events generated using Herwig and Jetset, respectively. The results for the 
gg and ginci. samples are seen to agree well over the entire jet energy range from 
5 GeV to 5 000 GeV, i.e. (n c h.)g incl ~ | ( ? ±/i.)ggevents for both Herwig and Jetset. This 
establishes that the properties of the g inc i. sample derived from the e + e~ data do 
indeed correspond closely to those oi gg events generated from a color singlet point 
source, confirming the results of ref. ta. 

Below, the results of a recent study £3 which uses this analysis method are pre¬ 
sented. The data sample employed is that qL the OPAL detector at LEP. Details 
of the event selection are presented in ref. £3. In total, 278 events are selected 
for the final gluon jet g; nc i. event sample. With the final cuts, Jetset with detec¬ 
tor simulation and the same selection criteria as are applied to the data predicts 
that both jets in the tagged hemisphere are quark jets with (83.0 ± 1.7)% proba¬ 
bility, where the uncertainty is statistical: this is the estimated purity of the g inc i. 
gluon jet sample. The mean energy of the gluon jets, (E , ) gincl , is determined to be 
(E )gind = 39.2 ± 0.3 (stat.) ± 1.8 (syst.) GeV. The uds quark jets in the study are 
defined inclusively using particles observed in event hemispheres opposite to those 
containing identified uds jets. The selection criteria are described in ref. £3. In total, 

28 007 uds hemispheres are tagged. The estimated uds purity of the sample, obtained 
by treating Jetset events with detector simulation in the same manner as the data, is 
(93.2 ± 0.2)%, where the uncertainty is statistical. The energy of the uds jets is given 
by the beam energy, 45.6 GeV. 

2 Results 

The mean charged particle multiplicity value measured for the gluon jets is (n c h .) gincl = 
14.63±0.38 (stat.)±0.59 (syst.). The corresponding result for the uds jets is (n c /j.) u dshemis. 
10.05 ±0.04 (stat.) ±0.23 (syst.). Before forming the ratio between the gluon and uds 
jet measurements, it is necessary to account for the different energies of the two sam¬ 
ples: the gluon jets have a mean energy of 39.2 GeV while the uds jets have a mean 
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Figure 1: (a) Predictions of QCD Monte Carlos for the mean charged particle multiplicity in gg 

event hemispheres and in e + e _ ginci. jets, as a function of jet energy. The corresponding predictions 
for uds qq event hemispheres are also shown, (b) Measured result for the multiplicity ratio r c h. 
between gluon and uds jets at -Ej et =39 GeV in comparison to analytic and Monte Carlo predictions 

shown as a function of -Ej e t- 


energy of 45.6 GeV. To account for the larger energy of the uds jets, the QCD analytic 
formula for the evolution of the mean event multiplicity in e + e _ annihilations Ed is 
employed. Assuming n/=5, the QCD evolution formula predicts the mean multiplic¬ 
ity in 78.4 GeV events to be (6.2 ± 0.4)% smaller than in 91.2 GeV events, where 
the uncertainty results from the maximum variation found by using the jet energies 
(39.2 and 45.6 GeV) rather than the event energies, n/=3 rather than rif— 5, and 
varying the value of as within its allowed range El Applying a multiplicative correc¬ 
tion of 0.938 ± 0.004 to the 45.6 GeV uds jet measurement presented above yields 
(rich.) =9.43 ± 0.06 (stat.) ± 0.22 (syst.) for the mean multiplicity of 39.2 GeV 
uds jets. 

The result for the multiplicity ratio r c h. between 39 GeV gluon and quark jets is 
therefore found to be: 

r ch. = 7 ^ C ' fe 39 , 2 G C eV = 1 ' 552 ± °- 041 ( stat 0 ± 0.061 (syst.) . (1) 

\ n ch.) u ds hemis. 

The systematic uncertainty is discussed in ref. 0. ft is to be noted that one of the 
systematic variations consists in using the Cone or JADE-E0 jet hnders to tag quark 
jets for the gi nc i. identification, rather than the k± jet finder: only a small change 
of Ar c h^0.02 occurs as a result of this variation. Thus, the result for r c h. changes 
by only about 4% if different jet hnders are employed for the analysis, in contrast 
to the Y events for which a difference of about 60% is observed (Table [1). This 
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emphasizes that the result ([[]) is only weakly dependent on a jet finding algorithm. 
This measurement of r c h. is shown by the solid point in Fig. [l](b). The experimental 
statistical uncertainty is indicated by the small horizontal bars. 

Various analytic results exist for the ratio r = (n)gi uon /(u) qua rk of the mean number 
of partons in a gluon jet to that in a quark jet. The original results B, valid to 
leading order, predict r to be r = Ca/Cf= 9/4. Later, higher order corrections 
valicLto the next-to-next-to-leading order were found to reduce this result by about 
1O%0. These results do not incorporate energy conservation into the quark and gluon 
branching processes. Recently, r has been calculated including not only the next-to- 
next-to-leading order terms but also energy conservations, and is found to be reduced 
yet further in magnitude. Momentum conservation is not included in this latter 
calculation, however: therefore energy-momentum conservation is only approximate. 
The analytic results for r, valid to leading order (l.o.), to next-to-next to leading 
order (n.n.l.o.), and including approximate energy-momentum conservation (n.n.l.o., 
E-cons.) are shown in Fig. |l](b) as a function of E- ]et = E cm J 2. For the evaluation 

of the strong coupling constant, as, the values nj =5 and A^hr o, =0.209 GeV0 have 
been used, where rif is the number of active quark flavors. The results for the two 
n.n.l.o. calculations are shown as bands: the upper edges of the bands show the 
results if the energy scale used to evaluate as is taken to be E cm ., while the lower 
edges show the results if this energy scale is taken to be E c m / 4. The widths of the 
bands therefore indicate the level of uncertainty associated with the ambiguity of the 
energy scale at which to evaluate as- The predictions of Herwig and Jetset at both 
the parton and hadron levels are included in Fig. |l](b). The width of the bands for the 
Monte Carlo results shows the variation which occurs if the shower cutoff parameters 
are allowed to vary within their allowed ranges 

For jet energies of 39 GeV, the n.n.l.o. calculation which incorporates energy 
conservation0 predicts values of r between 1.83 (if nf— 3 and the energy scale of as 
is E cm ) and 1.64 (if Uf= 5 and the energy scale is E cm / 4): this last value is only 
slightly above the measured result given above in relation (|l|). The analytic result is 
valid for quarks and gluons while the measurement refers to charged hadrons. Jet- 
set predicts a hadronization correction for r c h., defined by the ratio of the parton 
to hadron level curves in Fig. |l](b) (for E j e t = 39 GeV), of 0.91. The corresponding 
prediction from Herwig is 1.02. In this sense, the hadronization correction can be 
estimated to be about unity and to have an uncertainty of about 10%. Given the am¬ 
biguities of the energy scale at which to evaluate as, of the number of active flavors n/, 
of the hadronization correction, and due to the approximate nature by which energy- 
momentum conservation is included, is it seen that the analytic calculation of Drcmin, 
Hwa and NechitailoO is in general agreement with the measurement. In contrast, the 
analytic results which do not incorporate energy conservation 013 are seen to be in clear 
disagreement with this measurement, even considering the theoretical ambiguities. 

3 Summary 

The Erst quantitative test of analytic predictions for the ratio of the mean particle 
multiplicity between gluon and quark jets has been presented. The technique that 
allows this test is the identification of gluon jets in an inclusive manner, by all the 
particles observed in the eventnhemisphere opposite to a hemisphere containing a 
tagged quark and antiquark jet 0. The resulting definition is in close correspondence 
to the definition of gluon jets in analytic calculations, for the Erst time in the analysis 
of high energy e + e _ data. In this study, the gluon jet measurement, valid for 39 GeV 
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jets, is compared to the corresponding measurement from light quark (uds) jets, which 
has also been defined inclusively. 

The result for the ratio r c h., the mean charged particle multiplicity of gluon jets 
divided by the corresponding value for uds quark jets, is r ch = 1.552 ± 0.041 (stat.) 
±0.061 (syst.). This result is substantially smaller than the predictions of analytic 
calculations which do not include energy conservation in the parton branchings. A 
recent analytic calculation B which incorporates approximate energy-momentum con¬ 
servation predicts a parton level multiplicity difference between 39 GeV gluon and 
quark jets in the range from about 1.64 to 1.83, depending on the choice for the 
energy scale Q at which the strong coupling constant as(Q) is evaluated and on the 
number of active quark flavors, rif. This latter prediction is in overall agreement with 
the measurement, given the uncertainties due to the approximate nature of energy- 
momentum conservation in the calculation, missing higher order terms, the energy 
scale, the number of active flavors rif, and hadronization. 
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